Mixedness and entanglement for two-mode Gaussian states 
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We analytically exploit the two-mode Gaussian states nonunitary dynamics. We show that in 
the zero temperature limit, entanglement sudden death (ESD) will always occur for symmetric 
states (where initial single mode compression is zo) provided the two mode squeezing ro satisfies 
< ro < | log(cosh(2zo))- We also give the analytical expressions for the time of ESD. Finally, we 
show the relation between the single modes initial impurities and the initial entanglement, where 
we exhibit that the later is suppressed by the former. 
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I. INTRODUCTION 

The concept of entanglement , a typically quantum mechanical property, is a natural consequence of the superposition 
principle for composite systems. It was discussed by Schrodinger who immediately realized its (at the time) seemingly 
"unphysical" consequences [TJ. Nowadays the growing interest in the subject is related to quantum information 
theory and potential technological applications |3] • Thus a lot of effort has been put both in the quantification of 
entanglement and in studying the consequences of deleterious environmental effects on this property |4]. It is well 
known today, both theoretically and experimentally, that in general such effects tend to destroy quantum properties; 
in the case of systems with one degree of freedom this happens only asymptotically [5HZ]. The studies concerning 
the degradation of quantum effects are vast in the literature, specially those related to continuous variables systems 
subjected to noisy channels, see |8j[9]. Some recent works devote their attention to the robustness of Gaussian and 
non-Gaussian states under dissipative channels |10| . 

It was recently realized that in bipartite systems that entanglement may disappear suddenly, phenomenon called 
therefore "entanglement sudden death" (ESD) or finite-time disentanglement fT2"] . The phenomenon is strongly 
related to geometrical properties of the set of mixed quantum states |13j . but from a physical point of view, the issue 
is still far from being closed. 

In the present work we study more specifically the entanglement properties of two mode Gaussian states under a 
nonunitary evolution. In ref. [5HZ] the entropy growth of single-mode Gaussian and non-Gaussian states coupled to a 
reservoir has been presented in analytical form. We show that in the two-mode case, single-mode squeezing plays an 
important role in entanglement dynamics, even for fully symmetric channels (phenomenon that appears also in the 
case of qu-bits [14 ). For the specific case of a zero temperature reservoir and symmetric states (where initial single 
mode squeezing of both modes is zq) we are able wo show that ESD will always occur if < ro < \ log(cosh(2zo)), 
where r is the two mode squeezing. We also show that there exists an upper limit for the degree of mixedness of a 
state so that it can exhibit entanglement. This result may turn out useful for experimental realizations of two-mode 
entangled Gaussian states. Finally, for symmetric and pure states evolving in a reservoir at zero temperature, we 
analytically present the time of ESD. 

The paper is organized as follows: in section [TT] we present some well known properties concerning Gaussian States; 
in the next section we briefly review the analytical dynamics of two-mode Gaussian states; in section [TV] we present 
our results. Finally, some conclusions are drawn. 

II. GENERAL PROPERTIES OF TWO-MODE GAUSSIAN STATES 



A two-mode Gaussian state with vanishing averages (a*) = tr(aj/?) = 0, i = 1,2, can always be written as 

p G = Si(zi, Z2)S 2 (r)o(vi,i>2)Sl(r)Sl{zi, z 2 ), 
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where Si(zi, zq) is the single mode squeezing operator. S^iv) is the two-mode squeezing operator and <j{v\. v 2 ) is the 
two-mode thermal state. More explicitly we have: 



S 1 (z 1 ,z 2 ) = 
S 2 (r) = 

o(vi,v 2 ) = 



exp 



k - o 



exp 



r(a\a\ - aia 2 )l , 



exp 

1 1 

vi + 1 v 2 H 



fx + 1 



^2 
^2 + 1 



(2) 
(3) 



|n) (n| (g> |m) (m| 



where Zi is the single-mode squeezing parameter of the mode i, is the annihilation (creation) operator of the i— th 
mode, r is the two-mode squeezing and Vi is the "mixedness" of the i— th mode (by "mixedness" we mean the initial 
number of thermal photons of the state), related of the thermal two-mode state a. We assume (a^) = in this work, 
since the entanglement properties are independent of them, and in our equations of motion the second momenta 
evolution decouple form the first momenta. We also choose, without loss of generality in this case, the squeezing 
parameters Zi and r to be real. 

This state is entirely described by the parameters given above. However, in order to handle entanglement properties, 
it is most convenient to write these parameters in terms of the corresponding covariance matrix (CM): 
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where rii = (aja^), mi = — (of) 
of an observable £. The CM can also be written as 
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(aia 2 ), and ($) denotes the quantum expectation value tr(£,o) 
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where Vi is a 2 X 2 matrix related to the mode i, and C is a 2 x 2 matrix that gives the correlations (both quantum 
and classical) between the modes. For later use, we define some invariants of the covariance matrix as [71 115j: 

Iv = dctVp, 

h,2 = detVt,2, 

I 3 = detC, 

h = tr [V^CZVsZC^Z] . 



(6) 



These quantities are invariants under local unitary operations and Z = diag{l,— 1}. Next we give the explicit 
connection between the parameters of the Gaussian state and the matrix elements of the covariance matrix: 
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where we have defined 
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Once the evolution of the covariance matrix is obtained, the evolution of the parameters of the state can be inferred. 

Since we are working with Gaussian states, there are necessary and sufficient criterion to determine if the state is 
entangled |16fiTH] . Simon [TB] has shown that for any two-mode Gaussian state, if the following inequality is observed 

S {Vp) = hh + (1/4 - |/ 3 |) 2 -h- l/4(/i + h) > 0, (12) 
the state is separable. For the purposes of this work, Simon's criterion is enough to study the entanglement dynamics. 



III. NONUNITARY DYNAMICS AND ITS ANALYTICAL SOLUTION 



The degradation of the quantum information content of Gaussian states is a subject of interest, both for the 
technological and/or experimental applications as well as for fundamental quantum mechanics in what concerns the 
classical limit. 

The usual approach to non-unitary dynamics is by means of master equations, which has found several successful 
applications in other areas of physics. Our master equation reads 

P = Cp, (13) 

where £ is a superoperator given by: 

C = ^2 {li{ n f + l)(2ai •at — a\di • - • at a*) + 7 l nf (2aJ • a* - a<at • - • a,at)^ . (14) 

i=l,2 

In the equation above, 7$ is the dissipation constant of the reservoir related to the mode i, nf is related to the 
temperature of the thermal bath of the mode i, and at( a ») ^ s the creation (annihilation) operator of the respective 
mode. 



The equation (14 1 models a linear coupling between the state (the two-mode Gaussian state in our case) with a 
thermal bath of quantum harmonic oscillators. The evolution of each term of the covariance matrix, evolving under 
the dynamics described above, is given by: 

m = e" 2 *^ ((-1 + e 2t ^) nf + cosh(2z°) (i£ cosh 2 r + (1 + i/g) sinh 2 r ) + sinh 2 z°) ; (15) 



m, = — c 



(f° - "1 + C 1 + v a + v k) cosh(2r )) cosh z° sinh z°; (16) 
l e -t(7i+72)(i + v\ + vl) cosh(2? + 2§) sinh(2r ); (17) 
.I e -t(7i+72)( 1 + ^0 + ^0) sinh (2r ) sinh(2r? + z§). (18) 
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In the equations above, i — 1 or 2 (mode 1 or mode 2) and k ^ i. The parameters are such that z° is the initial 
single mode squeezing, 7^ is the dissipation constant, vf is related to the initial mixedness and nf is the reservoir 
temperature, where all the quantities refer to the i— th mode, as denoted by the suscript i. Also, we have that r$ is 
the initial two-mode squeezing. 



IV. RESULTS 



A. Entanglement dynamics and single mode squeezing 

In order to get a clear picture and to gain physical insight into the rich and multifaceted aspects of the non-unitary 
dynamics of general two-mode Gaussian states, we consider firstly the simplest case, i.e. the two-mode squeezed 
vacuum in dissipative channel. 

Let us consider the case of a two-mode vacuum state without single-mode squeezing (i.e. z% = z<z = 0). As shown in 
Figure [TJ for symmetric and asymmetric channels, with temperatures nf and nf different from zero, there will always 
be a finite time when entanglement vanishes. This can be understood using a geometrical picture of entanglement 
decay |13| . In fact, in this case, the long-time state is a separable mixed state, well within the set of separable states. 
Thus, if an initial state is entangled, it necessarily crosses the border of separable states in finite time. For zero- 
temperature baths, i.e. if nf = nf = 0, even if one uses different dissipation constants 71 and 72, the entanglement 
only disappears asymptotically. 
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Figura 1: Simon criterion for several cases of initial parameters. The common parameters are: z° = = v i = v 2 = 0, ro = 1. 
In each curve: gray: 71 = 0.1,72 = 0.5, m = 0.2, 122 = 0.2; blue: 71 = 0.1,72 = 0.1, ni = 0.2, ni — 0.2; green: 71 = 0.5,72 = 
0.5, m = 0.0, n 2 = 0.0; red: 71 = 0.1,72 = 0.1, m = 1.0, n 2 = 0.0; black: 71 = 0.1,72 = 0.1, m = 0.0, n 2 = 0.0; pink: 
71 = 0.1, 72 = 0.5, ni = 0.0, 712 = 0.0. 



Next we introduce single mode squeezing, i.e. z\ 7^ and/or z® 7^ 0. We note in Figure [2] that, even for the zero 
temperature case, one observes entanglement sudden death (ESD). We note that there is a dynamical entropy increase 
of the squeezed mode, caused by the reservoir, which acts in such a way that, for a relatively short time interval, this 
mode's entropy grows and then decays to the vacuum. We have observed that this dynamical entropy growth [SJ |B] 
causes the entanglement disappearance in finite time. 




Figura 2: Simon criterion for two set of initial parameters, say: v° — ^2 = 0,ro = 1. Blue: 71 = 72 = 0.1, nf = nf = 0, 2? = 
z% = 0; Red: 71 = 72 = 0.1, nf = nf = 0, z\ = % z% = 0. 

Since single-mode squeezing turns out to play a significant role on ESD. In Figure [3] we show Simon's criterion 
evolution for states with compression in only one of the modes, for reservoirs at null temperature. The vertical axis 
corresponds to the initial parameter z\ and the plot shows, for each instant of time, whether the evolved state is 
entangled (shaded area) or not (blank area). For instance, if the initial state has z\ = 2, Simon's criterion will change 
from negative to positive in finite time; if one have z\ = 1.0, the entanglement will vanish asymptotically. We will 
show an analytical relation between single-mode squeezing and ESD hereafter. 

In the following we want to discuss the instant of time in which the state becomes separable, or when occurs ESD, 
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Figura 3: Simon criterion for values of the initial parameter 2? and time t, where we have used z\ = z\. For the shaded area 
the Simon criterion S is negative (entangled state), in the white area S is positive. Note that there is a limit for z\ in which 
the state will have ESD. In this case: Z\ ~ 1.4. Parameters: ro = 1, v\ — u° — 0, nf = n§ = 0, 71 = 72 = 0.1. Graphs (a) and 
(b) are the same function, in a different time scale. 



in the case of zero temperature and symmetrical states (z® — 0, r > 0, v^—Vq and 7^ = 7, nf — 0, for i = 1,2). In 
this case the elements of the covariance matrix, which do not vanish, depend on 

n = n(t) = i e - 27 *((2^ + 1) cosh(2r ) - 1) = n e~ 27 *, (19) 

m = m(t) = ie" 27t (2i/ + 1) sinh(2r ) = m e _27t . (20) 

In this case Simon S can be factorized as S = (m + n){m + n + l)(m — n){m — n — 1). Since the first two factors 
are clearly positive, we can see that S is negative ]£ n < m < n + 1. If this inequality is satisfied at the initial 
time, no < m o < n o + 1- Now, multiplying by e~ 27 * we get n < m < n + e~ 27 * < n + 1. Thus, if the state is 
initially entangled, the evolved state is also entangled for any finite time. Since Simon S vanishes asymptotically 
lim^oo S(t) = 0, we see that either the entanglement decay is asymptotic, or the initial state is already separable. A 
two-mode squeezed vacuum, Uq = 0, is always entangled; hence, it separates asymptotically. 

Now we study tESD for states containing single-mode squeezing, i.e. z\ = z<i. Here one can find for tssD- 

e -2 7 t ESD _ 2 exp(r ) cosh 2z sinh r - 2 sinh 2 z 
exp(2ro)(cosh2ro — sinh2zo) 

where we consider r > 0, and z\ — z\ = Zq- In terms of the new variables r\ = exp(2r ) and £ = exp(2z ), the 
disentanglement time 

-2 7fESD _ 77(1 + C 2 - 20 
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is much easier to analyze. This equation has a valid solution when the right-hand side varies between and 1, that 
is, when 77 satisfies the inequality 1 < rj < | \ C + j^j- Going back to the original variables we conclude that the 

initial state separates at a finite time if < ro < \ log (cosh(2zo)). In this case, one can see clearly that the effect of 
single-mode squeezing is crucial to determine when the entangle will vanish. 



B. Effects of the mixedness in the initial state 



Recently |5J |5] it has been shown that in the case of single mode Gaussian states, there is an upper limit for the 
degree of global purity (represented by uf) of the state above which no quantum properties are visible. We now show 
that an analogous result holds for two mode Gaussian systems also. Following the steps in ref. O E] we can show 
that the initial state is entangled if the following inequality is satisfied: 

1 cosh _ x HI + A? ± 24(1 ± y ±4*1) + K) 2 (l + 8,0(1 + ,Q))X 
4 V (l + ,° + ^ ) 2 J 
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In Figure [4] we show the entanglement in the initial state as measured by S, with ro = 1. Notice that there is an 
upper limit on the initial state impurity above which the state becomes separable. 




Figura 4: Simon criterion for the two-mode Gaussian state, equation M. Note that impurities values sufficiently high can 
suppress the effect of entanglement given by the two-mode squeezing parameter r. Parameter: ro = f . 



V. CONCLUSION 



In this paper we review some aspects of two-mode Gaussian states, showing both the evolution of the state parame- 
ters (that entirely characterize the state) under nonunitary evolution and how entanglement of the state, characterized 
by the Simon criterion (12 1, depends on the initial state parameters. We show that, even in completely symmetrical 
reservoirs and zero temperature, the entanglement can vanish in finite time, depending on the single-mode squee- 
zing of the state. We give a condition for ESD relating single mode squeezing and two mode squeezing, that is 
< ro < \ log(cosh(2zo)). We analytically present the time when occurs ESD for symmetrical states, evolving in a 
reservoir with zero temperature. We also show that entanglement can be "suppressed" by the initial mixedness of the 
modes, iq: the initial two mode squeezing has an upper limit as a function i/f, and for values above this limit the 
state is separable. This can be helpful for experimental procedures, since any state will have some minimal impurity. 
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